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Abstract 

' Theoretical predictions for the production cross sections of leptoquarks at the CERN LHC are 

presented including higher-order QCD corrections. These corrections reduce the dependence on 
^ , the renormalization/factorization scales significantly. Moreover, they are required to exploit the 

H ' leptoquark mass reach of the LHC experiments. In this sequel to an earlier analysis performed for 

the Tevatron collider we extend the leptoquark analysis to the LHC energy. 
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Leptoquarks [l| have been searched for in the past at all high energy colliders operating 
at the energy frontier. If the Yukawa couplings for scalar leptoquark couplings to quarks 
and leptons of the first and second generation [on which we will focus in this report] is 
kept at small values, as required by low-energy precision experiments 0], the most stringent 
model-independent bounds on leptoquark masses, for branching ratios BR = 1 to charged 
leptons, have been set with M^q > 230 GeV and 241 GeV by the FNAL Tevatron CDF 
experiment for first and second generation leptoquarks, respectively, and correspondingly 
with 238 and 186 GeV by the DO experiment as reported recently in Ref. j^; these recent 
individual values come close to the earlier combined limit of 242 GeV for the first generation 
in Ref. j^. The DO and CDF experiments have also searched for leptoquarks decaying into 
quark and neutrino final states, resulting in a lower mass limit of 95 GeV and 117 GeV, 
respectively, in these decay channels Bounds of 290 GeV for the first generation could be 
set by the Zeus jsf and the HI experiments in direct electron-quark formation at DESY's 
HERA collider for Yukawa couplings of size 0.1 and electromagnetic strength, respectively. 
A limit _of 300 GeV is expected to be reached finally also by the Tevatron at the end of 
Run II 



The search for these novel particles will be continued soon at the CERN LHC. Preliminary 
feasibility studies by the LHC experiments ATLAS [9] and CMS [lo| indicate that clear 
signals can be established for masses up to about Miq ^ 1.3 to 1.4 TeV for first and second 
generation scalar leptoquarks, with a final reach of presumably 1.5 TeV. 

In this Brief Report we present the cross section for pair production of scalar leptoquarks 
at the LHC including next-to- leading order QCD corrections.^ The results are based on the 
calculation presented earlier in Ref. jlll] for Tevatron energies. The higher-order corrections 
must be included in the theoretical predictions for the production cross sections in order 
to reduce the strong dependence on the renormalization and factorization scales which the 
lowest-order Born term calculation is bedevilled with. Moreover, for the standard choice of 
the common renormalization/factorization scale near the leptoquark mass, the cross section 
is increased by the QCD corrections and the mass range that can be probed experimentally, 
is extended correspondingly. 

The basic processes for the production of leptoquark pairs at the LHC are gluon-gluon 
fusion and quark-antiquark annihilation: 

g + g ^ IQ + IQ 

q + q^IQ + IQ (1) 

as shown in Fig. ^ In the pointlike limit, the gluon-leptoquark interactions are determined 
by the non-abelian SU(3)c gauge symmetry of scalar QCD so that the theoretical predictions 
for the pair production of scalar leptoquarks are parameter-free in the first two generations 
1 terms can safely 

C =^\^ (41 

a]:^=^lp\ (2) 



for which Yukawa terms can safely be neglected 



para 
0: 



p (41 - 31/3') + (I8/52 -(3^- 17) log - 



P 



Note that QCD corrections to the production of vector leptoquarks are not under proper theoretical 
control as vector leptoquarks are described by a non-renormalizable effective theory. 
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where = (1 — 4M|q/s)^/^ denotes the leptoquark velocity and s^/^ the invariant energy of 
the subprocess. As expected, for low leptoquark masses gluon fusion is by far the dominant 
production mechanism. Quark- ant iquark annihilation becomes important for larger masses, 
providing 30% of the production cross section for leptoquark masses of about 1.5 TeV, as 
will be proved later in detail. 

The calculation of the QCD radiative corrections for the total cross section of leptoquark 
production 

a[pp LQ + LQ\ = agg + (Tgq + (Tgq (3) 

including virtual effects, gluon bremsstrahlung and initial-state parton splittings at 0(a;s) 
has been described in Ref. [l3| in detail. In the following we define the calculational 
scheme and the physical input parameters chosen for the present numerical analysis. The 
renormalization of the strong coupling as{^Ji) and the factorization of initial-state collinear 
singularities at /i are performed in the MS scheme. The top quark and the leptoquark are 
decoupled from the running of as(/^)- For the calculation of the pp cross section we have 
adopted the CTEQ6L1 and CTEQ6M jl^] parton distribution functions at LO and NLO, 
corresponding to = 165 MeV and A^^ = 226 MeV at the one- and two-loop level of the 
strong coupling a;s(/i), respectively. 

The QCD corrections strongly affect the parton cross sections near the production thresh- 
old, in particular in gg fusion which rises steeply with energy, (i) Sommerfeld rescatterin^ 
i.e. Coulombic gluon exchange between the final-state partons, is singular near threshold [15[ 
in the velocity of the produced leptoquarks, thus removing the phase-space velocity factor 
j3 from the cross sections for S-waves and damping the phase-space suppression — >■ 
for P-waves. These corrections are positive in the attractive color-neutral singlet channels 
and negative in the repulsive color-octet channels. For gluon initial states the color-singlet 
channels dominate, for quark-antiquark initial states the color-octet channels. The gluon- 
exchange ladders can be resummed giving rise to the threshold correction factors 
F// = 1 ± ixi,8 ^ - e-^i) and a;8/(e^» - 1) with xi = |^ and xs = |^ for singlet 

and octet channels, respectively, (ii) In addition, large initial-state gluonic bremsstrahlung 
corrections of the type log^/5 and log/3 emerge, where the universal double logarithm can 
be exponentiated. The perturbative expansion of the total parton cross section can be 
expressed in terms of scaling functions. 



where we define the Born term cross sections of the parton subprocesses by and the QCD 
virtual+soft and hard gluon corrections by f^+^/^ , The scaling functions depend on the 
invariant parton energy s^/"^ through (3. The expressions and /^+^/^//^ near threshold, 



1 99 3 h9 54 



fV+S/fB^ 11 fV+S/fB^ L (5) 

•'99 I J 99 3 36/3 "^"^ 48/5 
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can be derived from Ref. [13| by choosing the mass Miq as renormahzation and factorization 
scale. The scaling functions for qg initial states are neither linearly nor logarithmically 
singular at threshold. 

For high energies the NLO parton cross sections approach non-zero limits asymptotically 
in contrast to the ~ 1/s scaling behavior of the Born cross sections. The asymptotic NLO 
parton cross sections read in the nomenclatura of the preceding paragraph: 

H ^ 2159 ^ ^ 2159 

432007r ' ^'^^ 1944007r ' ^ ' 

The ratio between the gluon-fusion and Compton term is large, 9:2, which is a consequence 
of the large color charges of the gluons compared to the quarks. 

The scale dependence of the theoretical prediction is reduced significantly when higher 
order QCD corrections are included. This is demonstrated in Fig. |21 where we compare, 
for an intermediate leptoquark mass of 1 TeV, the dependence on the common renormal- 
ization/factorization scale /i at the leading and next-to- leading order of the total cross sec- 
tion. For the comparison of the LO and NLO results, we have calculated all quantities 
[i.e. the partonic cross sections, Q;s(/i) and the parton densities] consistently in leading and 
next-to-leading order. The scale dependence of the leading-order cross section is steep and 
monotonic. At next-to-leading order the scale dependence is strongly reduced, and the NLO 
cross section runs through a broad maximum near /i k, Mjjq/A which stabilizes the NLO 
prediction very effectively. 

The QCD radiative corrections enhance the cross section for the production of lepto- 
quarks in the vicinity of the standard scale choice fi ~ M^q, see Fig. |21 The i^-factors, 
K = ctnlo/clO) with all quantities in the numerator and denominator calculated in NLO 
and LO, respectively, are displayed in Table HI [Note that agg is only the remnant of the 
qg cross section after the collinear singularities are subtracted via mass factorization; over- 
compensation can give rise to negative values and agg must not be interpreted as a physical 
cross section.] The i^'-factors range from K ^ 1.5 at Miq ^ 200 GeV up to ~ 1.9 
at the upper end of the LQ mass reach of ~ 1500 GeV at LHC The leptoquark results 
coincide with the cross sections for the production of squark-antisquark pairs in the limit 
of large gluino mass 0, 0, 3] after the small (< 0.5%) contributions from rescattering 



diagrams involving four-squark self-couplings are eliminated. This comparison provides us 
with an independent check of the present numerical analysis. [Agreement with the scalar 
cross sections for low-mass leptoquarks in the compilation of Ref. [20] is only of limited value 
as the if-factors presented in Ref. ^] have been obtained by using the NLO code of the 
calculation Ref. [11].] Despite the sizable corrections at the scale /i ~ Miq, the moderate 
scale dependence of the NLO result implies a reliable perturbative expansion.^ 

To study the uncertainty in the prediction of the cross section inferred by the parametriza- 
tion of the parton densities, we have repeated the calculation using the MRTS2002 
parametrization ||23|]. The difference between the NLO results based on MRST2002 and 



^ For the production of Higgs bosons in gluon-gluon fusion, NLO X-factors have been found in the range 
between 1.5 and 2, cf. Ref. 



21 1, but they are modified only by another 20% in NNLO From this 

example we may infer that also the isomorphic perturbative QCD expansion in leptoquark production is 
expected to converge sufficiently fast and that the NLO results arc of physical significance. 
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CTEQ6, displayed in Table HI increases with increasing leptoquark mass due to the uncer- 
tainty in the gluon distribution at large x, but does not exceed 10% in the experimentally 
accessible mass range up to M^g ~ 1.5 TeV. 

In conclusion. As anticipated from many calculations performed in the past, the NLO 
analysis stabilizes the theoretical prediction of the production cross section for leptoquarks. 
Moreover, for the standard choice of the renormalization and factorization scales near the 
leptoquark mass, the higher order corrections increase the cross section so that Born cal- 
culations, performed for the same scale of QCD coupling and parton densities, provide us 
with a conservative lower limit. By the same token, the NLO corrections shift the LQ mass 
limit upward by an amount of about 100 GeV at the upper end of the LQ mass spectrum 
that can be probed at LHC 
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FIG. 1: Generic set of leading-order Feynman diagrams for leptoquark pair production through 
gluon-gluon fusion and quark-antiquark annihilation. 
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FIG. 2: Variation of the LO and NLO cross sections for pp LQ + LQ + X at the LHC with the 
renormahzation and factorization scales. The leptoquark mass has been set to Miq = 1 TeV. 




FIG. 3: Total cross section for pp LQ + IQ + X aX the LHC energy ^/s = 14 TeV as a 
function of the leptoquark mass Miq. The variation of the NLO cross section with the value of 
the renormalization/factorization scale is indicated by the shaded band. 
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TABLE I: Total cross section and i^-factors for pp LQ+LQ+X at the LHC energy \/5 = 14 TeV. 
The renormalization and factorization scales have been set to /i = Miq. The CTEQ6L1 and 
CTEQ6M IJ] parton densities have been adopted as default, whereas in the rightmost column 
a comparison is performed with the MRST2002 [2^ parton densities. The relative weight of gg, 
qq and gq contributions to the cross section at NLO is given in the second-to-last column. [The 
negative sign of Cgg is a mere artifact of subtracting collinear initial-state singularities via mass 
factorization; by definition, this particular higher-order quantity must not be interpreted as a cross 
section.! 



Mlq [GeV] 


CTEQ6(L0/NL0) 


MRST2002 


^LO [fb] 


o-NLO [fb] 


gg ■ qq 


gq 


K 


O-NLO [fb] 


200 


0.500 X 10 


2 


0.742 X 10^ 


0.94 


0.05 


0.01 


1.48 


0.779 X 10^ 


400 


0.140 X 10 


1 


0.224 X 10^ 


0.91 


0.10 


-0.01 


1.60 


0.243 X 10^ 


600 


0.135 




0.225 


0.88 


0.15 


-0.03 


1.67 


0.245 


800 


0.219 X 10" 


-1 


0.378 X 10^1 


0.84 


0.19 


-0.03 


1.73 


0.406 X 10"^ 


1000 


0.471 X 10" 


-2 


0.836 X 10^2 


0.82 


0.22 


-0.04 


1.77 


0.879 X 10^2 


1200 


0.121 X 10" 


-2 


0.221 X 10^2 


0.81 


0.24 


-0.05 


1.83 


0.226 X 10^2 


1400 


0.349 X 10" 


-3 


0.655 X IQ-^ 


0.79 


0.26 


-0.05 


1.88 


0.650 X 10^2 


1600 


0.109 X 10" 


-3 


0.210 X IQ-^ 


0.78 


0.28 


-0.06 


1.93 


0.201 X 10-3 


1800 


0.357 X 10" 


-4 


0.713 X 10-"^ 


0.77 


0.29 


-0.06 


2.00 


0.656 X lO"'^ 


2000 


0.122 X 10" 


-4 


0.253 X 10"^ 


0.77 


0.30 


-0.07 


2.07 


0.222 X lO"'^ 
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